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ABSTRACT 

In this paper we explore the potential effects of DM annihilations on the cosmological recom- 
bination spectrum. With this example we want to demonstrate that the cosmological recombi- 
nation spectrum in principle is sensitive to details related to possible extra energy release dur- 
ing recombination. We restrict ourselves to DM models which produce a negligible primordial 
distortion of the CMB energy spectrum (usually characterized as p and y-type distortions). 
However, since during the epoch of cosmological recombination (z ~ 1000) a large fraction 
of the deposited energy can directly go into ionizations and excitations of neutral atoms, both 
the cosmological recombination spectrum and ionization history can still be affected signifi- 
cantly. We compute the modifications to the cosmological recombination spectrum using our 
multi-level H i and He i recombination code, showing that additional photons are created due 
to uncompensated loops of transitions which are induced by DM annihilations. As we illus- 
trate here, the results depend on the detailed branching of the deposited energy into heating, 
ionizations and excitations. This dependence in principle should allow us to shed light on the 
nature of the underlying annihilating DM model (or more generally speaking, the mechanism 
leading to energy injection) when measuring the cosmological recombination spectrum. How- 
ever, for current upper limits on the potential DM annihilation rate during recombination the 
cosmological recombination spectrum is only affected at the level of a few percent. Never- 
theless, we argue here that the cosmological recombination spectrum would provide another 
independent and very direct way of checking for the presence of sources of extra ionizing or 
exciting photons at high redshifts. This would open an new window to possible (non-standard) 
processes occurring before, during and between the three epochs of recombination. 

Key words: Cosmic Microwave Background: cosmological recombination, temperature 
anisotropies, cosmological recombination spectrum, spectral distortions 



1 INTRODUCTION 



The 
and 



anomalies in the 
positron s at 1 GeV 



Adriani et al 



20091). 



Fermi (|Abdo et al 
H. E. S. S. Collaboration: F. Aharonianl 



cosmic-ray spectra of electrons 
- 1 000 GeV se e n us ing Pamela 
^20091) . H.E.S.S. 
2009IK and Atic 
Chang et alj l2008h ould be interpreted as signatures from 
annih i lating dark matter (e.g. seelHooper et al]l2004l : Icholis et al .1 
120081 : | Arkani-Hamed et all 120091 . and references therein). This 
very intriguing possibility has recently motivated several indepen- 
dent groups to reconsider the effects of such dark matter (DM) 
models on the ionization history of our Universe and the cosmic 
microwave b ackground (CMB) temperatur e and polarization 
anisotropies dGalli et alj 120091 : Isiatver et al.1 120091 : ICirelli et al.1 
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120091 : iHuetsi et alj 120091: iKanzaki et alj l2009h using modifi ed 
versions of Recfast dSeager et al.lll999ll200d : IWong et alj2008h . 



It is clear that the ionization history is sensitive to the number 
of extra ionizations and excitations of neutra l hydrogen or helium 
atoms at z ~ 1000 (e.g. IPeebles et al1l200d) . These extra ioniza- 
tions or excitations can be mediated by particles (e.g. electrons or 
positrons) or photons, which, for exa mple, could be produced as a 
conse quence of DM annihilations ( e.g. 



Padmanabhan & Finkbeine. 



nen 



1 20051) or decaying particles (e.g. IChen & Kamionkowskil 1200 
Therefore, both DM annihilations or decays of long-lived unstable 
particles in principle are able to delay reco mbination, introducing 
change s to the Thomson visibility function dSunvaev & Zeldovichl 
Il970al) at z ~ 1100, which then affect the cosmic microwave back- 
ground (CMB) temperature and polarization anisotropies, changing 
the position of the acoustic peaks and their (relative) heights. This 
allows one to place interesting constraints on possible DM annihila- 
tions during recombination using current and future CMB data (e.g. 
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IZhang et al.|[20o3 : iGalli et alj|2009l : Isiatver et alj|2009». Similarly, 
energy release by decaying particles l Zhang et ail 20071) . or more 
general sources of additional ionizations or excitations of neutral 
atom s during recombination feean et all 120031 120071 : IGalli et ail 
120081) can be constraint. Since the presence of any such sources 
could compromise our ability to measure the spectral index of the 
primordial power spectrum and its running, it is very important 
to consider these possibilities carefully, along with other physical 
corrections to the modeling of cosmologic al recombination (see 
iFendt et alj|2009l : ISunvaev & Chluball2009l. for detailed overview 
on recently considered processes). 

However, the CMB temperature and polarization anisotropics 
are not the only (direct) signals that can tell us about (non- 
standard) processes occurring during cosmological recombination. 
It is well known that the recombination of hydrogen and he- 
lium in the Universe leads to the emission of several photons 
per b a ryon, modifying the CMB energy s p ectrum ( Zeldovic h et al.l 
19681 : [Peebles! 1 19681 : iDubrovichl 1 19751 : bubrovich & Stolvarovl 



19971) . Recently, detailed computations of the cosmological recom- 



bination spectrum w ere carr ied out (e. g. Rubino-Mart in et al ■l200d : 
IChluba & Sunvaevl l2006bl : iRubino-Martni et al.l 120081) . showing 
that the recombinations of hydrogen and helium lead to relatively 
narrow spectral features in the CMB energy spectrum. These fea- 
tures were created at redshift z ~ 1300 - 1400, ~ 2100 - 2400 and 
~ 6000, corresponding to the times of H 1, He i and He ii recombina- 
tion, and, due to redshifting, today should still be visible at mm, cm 
and dm wavelength. Observing these signatures from cosmological 
recombination may offer an independent way to determine some 
of the key cosmological parameters, such as the primordial helium 
abundance, the number density of baryons and the CMB mon opole 
temperature at recombination (e.g. lChluba & Sunvaevl2008bh . Fur- 
thermore, it will allow us to directly check our understanding of the 
recombination process an d possible non-standard aspects (e.g. see 
ISunvaev & Chlubal2009l for a n overview), for example, i n connec- 
tion with early energy release dChluba & Sunvaevll2009dl) . 

In this work we will demonstrate that the cosmological re- 
combination spectrum also is sensitive to the branching of energy 
released due to DM annihilations into ionizations and excitations. 
As we show here, it is not only important how much energy is de- 
posited in total, but also when. Depending on the underlying model 
for the annihilating DM these efficiencies will differ, so that observ- 
ing the cosmological recombination spectrum may offer another 
very direct way for constraining such models. 

In earlier considerations of the possible effects in connec- 
tion with energy release during recombinat ion, this branchin 
was either param e trized with single numbers dPeebles et alj|200~ 
iBean et all I2003L 2007) , or s imple approximations were used 
llChen & Kamionkowskil l2004l: IPadmanabhan & Finkbeinenl2005l : 
iMapelli et al.ll2006h . A detailed account for all possible aspects of 
the problem related to a computation of these efficiencies is beyond 
the scope of this paper. Based on earlier investigations, we will 
therefore restrict ourselves to some simple examples (see Sect. [2jl, 
which are mainly intended to show the principle dependencies of 
the recombination spectrum on DM annihilations. However, one 
can carry out similar computations in connection with decaying 
particles or other, more speculative sources of extra ionizations 
during recombination. Since in those cases the cosmological re- 
combination should also exhibit signatures of these non-standard 
processes, we hope that this work will provide further motivation 
towards refined studies in connection with the cosmological recom- 
bination spectrum and potential future experiments measuring it in 
detail. 



But how do additional ionizations or excitations during re- 
combination actually affect the cosmological recombination spec- 
trum? Every extra ionization of hydrogen liberates an electron 
and proton. At z ~ 1000 this process is reversed by a recom- 
bination of the proton with another free electron after a rather 
short time. Since the H 1 L yman continuum is completely blocked 
dChluba & Sunvaevl [2007h . the electron is captured into some ex- 
cited state (n > 2), emitting at least two photons in the subsequent 
cascade towards the ground-state. This increases the total emission 
of photons by hydrogen during recombination, since in every addi- 
tional loop of transitions several quanta can be produced. Similarly, 
extra excitations allow additional electrons to reach high levels, so 
that in total more recombination photons will be released, by both 
hydrogen and at z ~ 2200 also by helium. 

The physics of this pro blem is very similar to the effect of he- 
lium photons on hydrogen dChluba & Sunvaevil2009fi) , leading to 
additional feedback-induced emission during recombination. Here 
the source of extra ionizations and excitations is related to photons 
emitted by helium at z ~ 2200 and ~ 6000 in the normal recom- 
bination process. Another example is connected with the changes 
introduced to the cosmological recombination spectrum as a con- 
sequence of ea rly energy release, which pro duces an primordial y- 
type distortion JZeldovich & Sunvaevl 19CT) of the CMB in the pre- 
recombinational epoch, leading to uncompensated loops of atomic 
transitions that attempt to restore the CMB blackbody spectrum 
(Chluba & Sunyaev 2009d). In this case, the extra ionizations and 
excitations are caused by the excess of photons in the Wien tail of 
the distorted CMB. 

For the cases considered here the latter process is not impor- 
tant, since the required energy that goes into additional ionizations 
or excitations during recombination can be tiny in comparison with 
the energy density of the CMB. In such cases the CMB spectrum is 
not affected significantly by the additional energy release through- 
out the entire pre-recombinational epoch (see Appendix [A]f or some 
estimates), while the dynamics of cosmological recombination can 
still be strongly modified. This is because the number of hydrogen 
and helium nuclei is a factor of ~ 2 x 10 10 smaller than the number 
of CMB photons. This makes it rather easy to perturb the recombi- 
nation process, while at the same time the CMB energy spectrum 
itself remains practically unaltered. 



2 MODELING OF THE DIFFERENT PROCESSES 

In this section we give a brief summary of the required equations 
for our multi-level recombination code to take the effect of DM an- 
ni hilations into account. W e close ly follow the approach outlined 
bv lChen&Kamionkowskil f2004h in connection with energy in- 
jection from decaying particles, inclu ding recent modifications and 
updat e s related to DM a nnihilations 1 Padmanabhan & Finkbeinen 
l2005l : IGalli et alj 120091 : Isiatver et al] 120091 : iHuetsi et al.l l2009i) 
However, some of the details here are slightly different, and we 
also introduce the basis for further improvements of our multi-level 
recombination code in connection with this problem. 



2.1 Overall energy injection rate 

Envisioning some self-annihilating DM particle x an d its antipar- 
ticle x, the total rate of energy release per unit volume is given by 
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dE 
dT 



= 2M x c 2 {av) N X N X 

* 2.9 x 1(T 31 [l+z] 6 eVs"'cm- 



100 GeV 



Q x h 2 



0.13 



(o-v) 



3 x 10- 26 cm 3 /s 



(la) 



(lb) 



Here M x = M x is the mass of the DM particle and its antipar- 
ticle; (crv) is the thermally averaged product of the cross-section 
and relative velocity of the annihilating dark matter particles; and 
N x = N x = N , x [l + z] 3 is the number density of DM parti- 
cles and their antiparticles, with the present day value N 0x » 

1.4X10- cm- 

Depending on the DM model, (crv) in general is a function of 
redshift. In particular, it could also include the effect of Sommer- 
feld enhancement dur i ng the epoch of cos mological recombination 
(e.g. lGalli et alj|2009l ; lsiatver et alj|2009h . which can be important 
since at that time the relative velocities of the DM particles be- 
come small. One can incorporate these possibilities by replacing 
(crv) = S(z) (<tv) , where (ay) = const, however below we will 
restrict ourselves to cases with (crv) = 5 (z) (<rv) = const. 

Another a spect of the probl em is connected with the clus- 
tering of DM fauetsi et alj|2009l) . This effect becomes only im- 
portant at low redshift (z < 100) and leads to an enhancement 

(N x N x \ = B(z) (Nx) °f me avera g e squared DM number density, 
with clustering boost factor B(z) > 1 . This effect can become very 
pronounced, depending on the assumed halo concentration model 
and the lower mass cutoff for the halo mass function. However, 
here we are mainly interested at the CMB spectral distortions gen- 
erated at z > 200, where the clustering of matter is negligible. We 
therefore neglect this aspect of the problem here, using B(z) = 1 
throughout, so that Eq. QJ remains unaltered. 

However, it is very important that depending on the involved 
annihilation channels (e.g. photons, leptons, hadrons, neutrinos) 
only a fraction, f A , of the released energy will be deposited into the 
intergalactic medium (IGM), going into heating, and ionizations 
or excitations of atoms (i.e. hydrogen and helium). For example, 
energy released in form of neutrinos (at redshifts of interest to us 
here) will be carried away, so that usually one expects f& < 1 . Fur- 
thermore, because the transparency of the Universe to photons and 
the energy deposition efficiency of different particles (e.g. electrons 
and positrons) depend on the redshift of injection and the cosmo- 
logical model (e.g. densities and expansion rate), f A is a function of 
time and cosmology. 

To include this aspect of the problem into the computations 
we therefore write 



d£d 
At 



= Mz) 



dE 
d7 



/d(z)eoA , H[l+z] 3 eVs- 1 , (2) 



with the dimensionless parameter 

,2 W 



e = 1.5 x 10 



-24 



M x c 



100 GeV 



n x h 2 



0.13 



(crv) 



3 x 10- 26 cm 3 /s 



(3) 



Here we also used A'h = 1.9 x 10 7 cnT 3 [l + z] 3 for the number 
density of hydrogen nuclei in the Universe. 

A detailed computation for f A as a function of time and 
cosmology i s beyo nd the scope of this paper. However, recently 
ISlatveretalJfe009n computed the function for different models 
of annihilating DM for the concordance model. Their models give 
typical values for f A ranging from / d ~ 0.1 to f& ~ 1, where the 



largest values are reached for annihilation channels XX ~ * e+e at 
high (z > 10 3 ) redshifts. They also provided some simple fitting 
formulae for specific DM models, which we will use below. For 
the purpose of this paper this should be sufficient. 



2.2 Heating of the medium 

As mentioned above, only part of the energy that is deposited into 
the IGM will go into heating of the medium. The rest will lead 
to ionizations or excitations of hydrogen and helium atoms. If we 
denote the fraction of the deposited energy that goes into heating 
by gh(z), then we can write the additional term in the evolution 
equation of the temperature of the medium which is related to DM 
annihilations as 



dt 



_2_ gh(z) dEj 

h~ 3k /V H [l+/He+*e] dt 



(4) 



Here A'h is the total number of hydrogen nuclei; /n e ~ 8% is the 
number of helium nuclei relative to the number of hydrogen nuclei; 
and X c = A'c/A'h is the usual free electron fraction. 

It is clear that at high redshifts, well before the epoch of re- 
combination, practically all the deposited energy goes into heating 
of the medium, so that g n (z) ~ 1. Due to energy conservation, one 
also has g h (z) = g d (z) - g,on(z) - gex(z), where here g ion (z) is the 
fraction of the deposited energy that goes into ionizations of atoms, 
and gcx(z) the fraction that goes into excitations. Note that every 
ionization event also leads to partial heating of the medium, since 
the liberated electron (and nucleus) will usually have some (large) 
excess energy which then will be dissipated, e.g. in form of sec- 
ondary particles, which again can lead to the heating and ionization 
of the medium. In detailed computations of the efficiencies g ion and 
g cx this has to be accounted for. Below we will specify which ap- 
proximations for these functions we use in this work. 



2.3 Ionizations and excitations of atoms 

Knowing that a fraction gi on (z) of the energy deposited by DM anni- 
hilations is going into ionizations of both hydrogen an helium, one 
can write g i011 (z) = g" n (z) + g" c n (z), where gjj n (z) and g*(z) are the 
partial contributions of hydrogen and helium, respectively. Simi- 
larly, one has g cx (z) = g"(z)+gj x = (z) for the energy that goes into ex- 
citations. Since the chemical mixture of the medium can vary, it is 
useful to introduce the specific contributions g\(z) = g\(z) N tot /Ni, 
with i = H or i = He; and a = ion or i = ex. For hydrogen one has 
Mot/Mr = (1 + /He) and for helium N to t/N We = (1 + / H e)//ife. 

For the net ionization rate from the ground states of neutral 
hydrogen and helium related to DM annihilations one then finds 



dN»' 



dt 

dA'" E1 



1 



.ft) 



dt 



1+/He 

/He 
1 + /He 



E H ' 

ion 

g He (z) 



d£ d 
dt 

d£ d 
dt 



(5a) 
(5b) 



where E Hl = 13.6 eV and £ Hei = 24.6 eV are the ionization poten- 

lon ion t 

tials of hydrogen and helium, respectively. For Eq. ((5) it was as- 
sumed that He hi is not important. Furthermore, it is assumed that 
the energy which is consumed in each ionization is equal to the ion- 
ization energy, so that j^- ^ | _ gives the rate of ionization events 
per unit volume. By definition of gj on this should be possible. 

With the same arguments, for excitations of hydrogen and he- 
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hum from the ground state one can write 



At 



1 gfc) d£ d 



1+/ Hc Eg' df 

/He ggO) d£j 
1+/ H e Eg? 1 df 



(6a) 
(6b) 



where now E^ = 10.2eV and a 21.0eV are the transition 
energies to the second shell of hydrogen and helium, respectively. 
Here in particular it is assumed that states with principle quantum 
numbers n > 2 are not directly excited. For the purpose of this paper 
this approximation will do, however, this may lead to an underesti- 
mation of the total number of additional secondary low-frequency 
photons produced in our computations. This is because excitations 
to highly excited levels (n » 2) will directly allow some electrons 
to make transitions among excited states, which leads to emission 
at low frequencies. On the other hand, for excitations to the sec- 
ond shell the electron will more likely stay within lower levels, and 
hence produce emission at high frequencies. 



2.3. 1 Expressions for g| on (z) and g^(z) 

Detailed computations of the specific ionization and excitation 
fractions gj (z) and gj, x (z) require to follow the evolution of pri- 
mary and secondary, non-thermal elect rons and photons pro duced 
by the DM annihilation process (see iKanzaki et all 12009. for a 
recent study in this connection). Such computations have to in- 
clude several cooling (e.g. Compton and Coulomb cooling for elec- 
trons) and particle creation processes (e.g. pair production by pho- 
tons), and various aspects of radiative transfer (e.g. photon feed- 
back, escape of photons from the main resonances). This is far be- 
yond the scope of this paper , however, based on calculation s by 
IShull & van Steenberd dl985l) . IChen & Kamionkowskil d2004l) pro- 
posed 



^ HHl 0) = lexO) 



i -x n 



(7) 



where X p = N p /Nn is the free proton fraction. In the con- 
text of cosmological recombination with DM annihilations 
this a pproximation has already b een u s ed by seve r al au- 



thors i Padmanabhan & Finkbeinej 20051 ; Zhang et al 



Mapelii et al. l l2006r . iGalli et al.l 120091: ISlatver et al 



2006 



2009 



Huet sietal.ri2009l) . Furthermore, IPadmanabhan & Finkbeinen 



1 20051) also applied a similar expression for the specific ionization 



„ i-z H ,n 



and excitation fraction of helium, i.e. g*(z) = g^fOO 
where Zh e „ = Nh eu /Nhc is the fraction of singly ionized helium 
atoms relative to the total number of helium nuclei. In this 
approximation g h = [1 + 2X P + / He (l + 2Z Heii )]/3[1 + / He ], where 
both ionizations and excitations were included. 

These are very rough approximations, since many details of 
the computations are not represented or recoverable in this way. For 
simple order of magnitude computations this approach certainly is 
acceptable, but for more detailed calculations that aim at including 
model-dependencies, refinements become necessary. For example, 
without further details it is not easy to say what fraction of excita- 
tions are due to photons and which are due to electrons or collisions 
in more general. In the former case, one should introduce modifica- 
tions to gj! x (z) due to photon escape, which will strongly depend on 
the actual density of hydrogen atoms and the expansion rate of the 
medium, but not only on the ionization fraction. For conditions in 
our Univers e, in that case it will be poss i ble to neglect excitations 
(e.g. like in IPadmanabhan & Finkbeinej d2005l )). since the escape 



probabilities, P esc , in the main resonances of hydrogen and helium 
are extremely small, so that g[! x (z) ~ P sx gl x (z) ~ 0. On the other 
hand, if excitations mainly occur due to interactions with primary 
or secondary electrons no such modification is n ecessary. At low 
energies (E < 100 eV) this seems to be the case JShull|[l979l) . As 
we will see below, the results do depend on this assumption, and a 
more detailed co mputation will be necessary. 

In addition, IShull & van Steenber 2 (1985) assumed that the 
ionization fractions of hydrogen and helium are always equal. This 
assumption is not valid in the cosmological recombination prob- 
lem, where helium is completely recombined at z ~ 1700 - 1800 
while hydrogen is still fully ionized. It is not clear that this extreme 
case can be obtained from their results, since details in the radiative 
transfer would be very different, li kely leading to non-line a r scal- 
ings. Already the computations of IShull & van Steenberd dl985l) 
suggest that the scaling of g?(z) with ionization fraction from neu- 
tral to fully ionized media already appears to be faster than linear. 
In the case of hydrogen they provide a fit that is close to 



(8) 



Also, their results suggest that for nearly neutral media one has 
gf on ~ 2/5, g» ~ 1/2, j» ~ 3/5 and g» e ~ 3/10 (see their Table 2 
with / He = 10%) instead of g" n ~ ~g» ~ g* ~ g«= ~ 1/3. 

All this shows, that for more accurate computations within the 
cosmological context, relatively large differences to the above ap- 
proximations for g"(z) can still be expected. A more detailed com- 
putation, where we also plan to include doubly ionized helium, ex- 
citations of levels with n > 2, and detailed radiative transfer, will be 
left for some future work. However, below we will demonstrate that 
the cosmological recombination spectrum is sensitive to the form 
of gj"(z), while for cases in connections with DM annihilations the 
dependence on f A is much weaker. 



2.4 Modifications to the multi-level recombination code 

For the computations presented in this paper we use our multi- 
level h ydrogen and helium recombination code dChluba & Sunvaevl 
l2009bl) . It in principle allows to include the effect of photons feed- 
back and several other recently considered physical processes that 
modify the recombin ation history of the Universe at the percent - 
level (for overview see lFendt et al ,l2009l : ISunvaev & Chlubal2009l) . 
However, here we will not include most of these corrections, since 
the ambiguities introduced due to DM annihilations are much larger 
in any case. We only want to demonstrate the principle aspects of 
the problem and show that the cosmological recombination spec- 
trum is sensitive to DM annihilations. 

As a first step one should add the term given by Eq. {4]l to 
the normal evolution equation for the temperature of the medium. 
Due to the tight coupling of the photon and electron temperature 
by Compton scattering, the additional heating will not affect the 
results of our multi-level recombination code until low redshifts 
(z £ 200), where the energy exchange between electrons and pho- 
tons becomes inefficient. However, the continuous heating of the 
medium at high redshift in principle will lead to some (small) 
primordial p or y-type spectral dist ortion of the CMB well be- 
fore the epoch of recombination (e.g.[ Zeldovich & Sunvaevl l 19691 : 
ISunvaev & Zeldovich|[l970bl ; lillarionov & Svunvaevlll975al lbl). We 
do not take this modification of the background radiation into ac- 
count, but according to our estimates (see Sect. [A}, for the DM 
models under discussion here, their annihilation should never lead 
to any important primordial CMB distortion (e.g. the y-parameter 
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Figure 1. Effect of dark matter annihilation on the free electron fraction for 
different constant values of fae®, as labeled. The curves in the upper group 
show the relative di fference in the free electron fraction in comparison to 
the reference model jRub ino-Maitm et alj200^) . The shaded area indicates 
the region around the Thomson visibility function, which defines the last 
scattering surface. 



Figure 2. Effect of dark matter annihilation on the number density of free 
protons, Np, and Heii ions, N\\ Ell , for different constant values of f&eo, as 
labeled. The shaded area indicates the region around the Thomson visibility 
function, which defines the last scatterin g surface. For comparison we also 
show the result for our reference model I Rubino- Martin et alj2008h . 



was always smaller than ~ 10 10 - 10 8 ). Therefore the effects dis- 
cussed in our previous work on pre-recombinational energy release 
(Chluba & Sunvaevll2009ci) are negligible here. Note that with re- 
spect to those earlier computations, the main difference is that we 
now include direct ionizations and excitations by DM-induced par- 
ticles, which, as mentioned in the introduction, can be very effi- 
cient. 

In the multi-level recombination code, one should also add 
Eqs. $5$ and (f6]l to the ground state rate equations of hydrogen and 
helium, accordingly. The additional ionizations introduced by DM 
annihilations will liberate an electron and (nucleus), and depending 
on the epoch at which this ionization occurs, the ionization will be 
directly compensated by the recombination of another electron. In 
the pre-recombinational epochs one should also allow for di rect re- 
combinations to the ground state dChluba & Sunvaevl2009rJ) . while 
during the recombination epochs electrons will be captured to ex- 
cited states (with principle quantum numbers n > 2), potentially 
liberating several (low frequency) photons in the cascade towards 
the ground state. Also, Eq. l |6al > should be subtracted from the rate 
equation for the H i 2p state in order to ensure conservation of the 
electron number. Furthermore, we subtract Eq. l |6bt from the rate 
equation of the 2'P[-level, assuming that with n = 2 only this level 
is reached. Since the transition rate to the 2 3 Pi-level is ~ 10 7 times 
smaller, this should be sufficient. Similarly, we subtract Eqs. l[5j 
from the rate equation for the free electrons. 



3 HOW DOES THE COSMOLOGICAL 

RECOMBINATION SPECTRUM CHANGE DUE TO 
DARK MATTER ANNIHILATIONS? 

In this section we illustrate the dependencies of the cosmological 
recombination spectrum and ionization history on the energy injec- 
tion by annihilating DM. We start the discussion with cases that 
assume f A (z) e Q = const, and then go to more complicated models 
in the following. In particular we want to demonstrate that the cos- 
mological recombination spectrum is sensitive to differences in the 



DM annihilation model, in particular related to the branching of the 
deposited energy into heating, ionizations and excitations. 

3.1 Effect of DM annihilation on N e : constant / d (z) eo 

Figure Q] shows the effect of DM annihilation on the free electron 
fraction. We included 10 shells for hydrogen and helium into our 
computations and assumed that all g*(z) are given by Eq. Q. One 
can clearly see that DM annihilations have the strongest effect at 
low redshifts, leading to a delay of recombination (at z ~ 1000) and 
an increase in the residual free electron fraction at very low redshift 
(z ~ 200). The recombination of neutral helium (at z ~ 2000) is 
hardly changing, even in the most extreme cases considered here, 
implying that the net recombination rate for helium is not affected 
as strongly by ionizations due to DM annihilation. 

One reason for this behavior is that for a given energy de- 
position rate, dE^I dt, due to the difference in the ionization po- 
tentials there are about ~ 2 times fewer ionizing and photons per 
helium atom available than for hydrogen (see Eq. (O for con- 
firmation). Similarly, the effective DM-induced excitation rate is 
~ 2 times smaller. Another reason is that any small relative differ- 
ence AN^ C /N^ c in the number of free electrons from helium, due 
to its small abundance (~ 8% in comparison to hydrogen), will 
have a ~ 13 times smaller effect on the total ionization history 
N c = N^} + N^ c , which includes the electrons from hydrogen. In 
the early stages of helium recombination one therefore expects that 
the ionization history can only be affected by a comparable amount 
as during hydrogen recombination when increasing the DM anni- 
hilation rate ~ 20 - 30 times. Looking at Fig.[T] and comparing the 
curves for / d (z) e = 2 x 10~ 23 and f d (z) e = 5 x 10~ 22 at z ~ 1 100 
and z ~ 2000, seems to confirm this statement. 

However, towards the end of helium recombination the main 
reason for the rather small effect of DM annihilations on the free 
electron fraction is connected with the acceleration of helium re- 
combination caused by the absorption o f resonant Hei photons 
in the Lyman continuum of hydrogen JSwitzer & Hiratal 120081 : 
iKholupenko et al.l200llRubino-Martin et alj200"8h . This shifts the 
end of helium recombination from z ~ 1600 to z ~ 1750, because 
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the effective recombination rate of helium is increased many times 
by this processes. It is extremely hard to delay helium recombi- 
nation with DM annihilations, once this process is working well 
(z s; 1900). Here it is also important that the Heii ions interact 
with a bath of free electrons from hydrogen. Per He ii ion there are 
about 13 electrons available for recombinations, while per proton 
during hydrogen recombination there is only one. This number of 
electrons per He ii ion remains practically constant until the recom- 
bination of hydrogen begins. Therefore, without DM annihilations 
practically all helium atoms recombinfl leaving basically no free 
Heii ions at low redshifts (see Fig.[2j- 

Only at very late stages, when recombinations of helium are 
already slow, one again expects some modifications in the helium 
ionization history. This is because there DM annihilations can (par- 
tially) reionize helium atoms, without this process being (signifi- 
cantly) reversed by recombinations. In Fig.[2]one can see that at low 
redshifts the number of residual He ii ions indeed increases strongly 
when accounting for DM annihilations. Nevertheless, in all consid- 
ered cases no more than ~ 0.3% of helium atoms are reionized 
by DM annihilation at low redshifts, and the total contribution of 
free electrons from helium to the residual free electron fraction at 
z ~ 200 does not exceed a few percent. 

3.2 Effect of DM annihilation on the cosmological 
recombination spectrum: constant / d (z) eo 

As we have seen in the previous section, the largest modifications 
in the free electron fraction appear at the end of hydrogen recombi- 
nation. For the Thomson visibility function and the CMB power 
spectra the modifications around z ~ 1100 are most important, 
while the huge relative changes in the residual electron fraction at 
z S; 500 actually do not matter that much. 

It is known that the recombination lines from hydrogen mainly 
appear at z ~ 1300 - 1400 (e.g. see iRubino-Martrn et al.l 120061 : 
IChluba & Sunvaevl2006al ; IChluba et al.l20071) . where about ~ 20% 
of the hydrogen atoms recombined, while at maximum visibility 
(z ~ 1100) already ~ 86% of all Hi was formed. From the differ- 
ences in the free electron fraction it is therefore already clear that 
the Hi recombination lines will mainly be modified on the blue 
sides of the recombination features, with an increase of the emis- 
sion due to additional ionizations and subsequent recombinations. 
With increasing DM annihilation efficiency the changes will be- 
come more strong and should eventually also affect the maxima 
of the recombination features. This will lead to shifts in their po- 
sitions towards higher frequencies, and an increase in the overall 
amplitude and width of the recombination features. 

In Fig. [3] we illustrate this behavior of the different compo- 
nent in the Hi recombination spectrum when taking the effect of 
DM annihilations into account. We included 10 shells for hydrogen 
and helium into our computations and assumed that all g-'(z) are 
given by Eq. 0. For all components the DM annihilation increases 
the overall amplitude of the emission and hence the total number 
of photons released during recombination. This is simply related 
to the fact that every extra ionization caused by DM annihilation 
will liberate an electron which then can recombine to some excited 
state of hydrogen. From there it will cascade towards lower levels 
emitting several photons on its way. These photons are released in 
addition to those from the normal recombination epoch. Because 

1 This would also be true without the inclusion of the speed up due to the 
H i continuum opacity. 
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Figure 3. Effect of dark matter annihilation on the H i bound-bound dipole 
recombination spectrum (upper panel), and free-bound and 2s- Is two- 
photon continuum emission (lower panel) for different constant values of 
/dfo, as labeled. Shown are the present day (z = 0) spectral distortion of the 
CMB generated by hydrogen during cosmological recombination including 
10 shells in the com putations. For comparison w e also show the result for 
our reference model jRub ino-Martm et al.l2008h . 

that at z < 2000 the Hi Lym an continuum is completely blocked 
(e.g. lChluba & Sunvaevll2007r) , every ionization caused by DM an- 
nihilations will lead to at least two photons (i.e. Balmer continuum 
and Lyman a) at lower frequency. However, since electrons also 
can be captured into some highly excited state (n > 2), the effec- 
tive number of emitted photons per ionization can be larger than 
two. In Fig. [3] one can see that indeed also the emission in transi- 
tion among highly excited levels, appearing at low frequencies in 
the cosmological recombination spectrum, increases when includ- 
ing the effect of DM annihilations. This shows that a significant 
number of electrons are captured to states with n > 2. 

Furthermore, from Fig.[3]it is clear that in particular for the Hi 
bound-bound dipole emission lines the positions and widths of the 
recombination features are affected by the DM annihilations. Both 
the free-bound and 2s- Is two-photon continuum emission are less 
sensitive in this respect, since they are initially very broad. In this 
context, especially the spectral features due to Paschen a (visible at 
v ~ 120 GHz) and Balmer a (visible at v ~ 350 GHz) are interest- 
ing, as they are both very prominent and not overlapping so much 
with other spectral features. Nevertheless, the distortions from the 
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Figure 4. Component to the total H i cosmological recombination spectrum 
which is induced by dark matter annihilation for different constant values 
of /dfo. as labeled. We computed the absolute difference, A/dm(v), in the 
total, p resent-day H i recombinatio n spectrum with respect to our reference 
model ( Rubino-Martfn et al. 21 108 1 which does not include DM annihila- 
tions (shown as thin solid line). 
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Figure 5. Effect of dark matter annihilation on the He i bound-bound dipole 
recombination spectrum for different constant values of /dto, as labeled. 
Shown are the present day (z = 0) spectral distortion of the CMB generated 
by neutral helium during cosmological recombination including 10 shells in 
the co mputations. For comparison we also show the result for our reference 
model iRubino-Martin et all2008l) . 



recombination epoch are affected in basically all spectral bands. 
From an observational point of view it will be important to look at 
the CMB distortions in many frequency channels and to determine 
the positions and width of several features simultaneously. 

Here we would also like to mention that the very high fre- 
quency spectral feature visible at v ~ 4 THz is created by the H i Ly- 
man a resonance at z ~ 600. Although from an observational point 
of view this distortion is not very interesting (the cosmic infrared 
background is far too strong in this spect ral band), for computat ions 
of the low redshift chemistry (e.g. see ISchleicher et al]|2008l for 
recent computations), a lso including the effect of non-equilibrium 
backg round radiation l lSwitzer & Hiratal 120051 : IVonlanthen et al .1 
l2009h . such feature may be relevant. However, in this case pro- 
cesses directly related to secondary particles from DM annihilation 
may still be m ore important (e.g. like for example found in the case 
of cosmic ravs lJasche et al .120071) . 

3.2.1 The Hi spectral distortion produced by DM annihilations 

For Fig. [4] we also separated the contribution, A/d m (v), to the to- 
tal H i cosmological recombination spectrum (bound-bound + free- 
bound + 2s- Is two-photon continuum) which is induced by dark 
matter annihilations. The main broad features at v ~ 440 GHz, 
150 GHz and 70 GHz in A/ dm (v) are due to reactivation of the 
Balmer a, Paschen a and Brackett a resonance at z < 1300 - 1400, 
respectively. The width of the associated spectral features is compa- 
rable to those produced in the normal recombination process. Fur- 
thermore, one can see that A/d m (v) oc f A e , so that it is possible to 
compute the total H i recombination spectrum by adding this com- 
ponent (with appropriate rescaling) to the normal recombination 
spectrum. For f A e ~ 10~ 22 the relative contribution of the DM- 
induced emission on average reaches a few percent, and in some 
bands up to ~ 10%. In the most extreme case presented here the 
relative change in the recombination spectrum introduced by DM 
annihilations reaches ~ 60% - 70% in the vicinity of the Balmer 
and Paschen series. 

Why does A/d m (v) actually exhibit spectral features? For this 
it is important that although DM particles continuously annihilate, 



the fraction of energy (z) °c [1 — X v ] that actually goes into ion- 
izations is very small until the recombination epoch of hydrogen 
is entered. Similarly, excitations of neutral atoms caused by DM 
annihilations are not effective until hydrogen starts recombining. 
On the other hand, the energy deposition rate per hydrogen atom 
scales like oc /d(z)[l +z] 3 , and hence decreases with redshift. 

Therefore, the product /d(z)[l +z] 3 g^Jx) exhibits a maximum, that 
is mainly determined by g? n (z) which changes rather fast close to 
the recombination epoch. As a consequence, the spectral distortions 
which are introduced by DM annihilations also show maxima with 
a width which is comparable to the duration of recombination. As 
we will see below (Sect. 13.4b . the DM-induced changes to the H i 
recombination spectrum are indeed sensitive to the redshift depen- 
dence of g? n (z) and to a smaller extend of / d (z). 

3.2.2 Effect on the helium recombination spectrum 

In Fig. [5] we also present the changes in the CMB spectral distor- 
tions introduced by Hei. Here the effect of DM annihilations is 
significantly smaller than for hydrogen: for f^co = 5 x 10~ 22 the 
differences are of the order of ~ 10% - 20% in some bands, while 
for hydrogen changes up to ~ 60% - 70% were found (cf. upper 
panel in Fig. [3}. Still this is much larger than the changes seen in 
the helium ionization history, where for / d 6 = 5 x 10~ 22 the correc- 
tions were ANf"/Nf e ~ 0.1% x 13 ~ 1% (cf. Fig.rjj. This implies 
that the recombination spectrum is more sensitive to energy deposi- 
tion during helium recombination than the ionization history itself. 
Also it is clear that the small changes in N e during helium recombi- 
nation will not propagate very much to the CMB power spectra, so 
that one cannot expect to see any signature of DM annihilation dur- 
ing helium recombination in the C/'s. However, directly observing 
the helium recombination spectrum in principle could shed light on 
processes occurring during this epoch. 

How does this work? In the early stages of helium recombi- 
nation and in its pre-recombinational epoch (z > 2600) the number 
of extra He i ionizations caused by DM annihilations is very small 
(g?5(z) oc [1 - Z Heii ] <k 1). According to our parametrization most 
of the deposited energy is going into heating of the medium at that 
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time (Sect. l2.3l ). The same is true for the pre-recombinational epoch 
of hydrogen, explaining why practically no extra emission is pro- 
duced, even though DM is continuously annihilating (Sect, [3~2Tl l. 

Our discussion in Sect. 13.11 has already shown that the ion- 
ization history of helium is not affected as much by DM annihila- 
tions, so that from this very few extra photons are expected. There 
it was most important that the small fraction of neutral hydrogen 
present at the end of helium recombination leads to a huge increase 
in the photons escape probability, so that helium recombination is 
strongly accelerated. This makes it very hard to change the He i ion- 
ization fraction by DM annihilations, since practically every extra 
ionization is directly reversed by a recombination. 

However, during the whole epoch of helium recombination, 
DM annihilations do lead to some extra ionizations. This drives 
loops of transitions, that start with the ionization of a neutral he- 
lium atom by DM annihilation, and end with the release of extra 
photons at lower frequencies in the cascade of electrons from ex- 
cited levels towards the ground state. This explains why it is pos- 
sible to see the effect of extra ionizations in the helium recombi- 
nation spectrum, while at the same time the effect on the ioniza- 
tion history is much smaller. However, due to the small abundance 
of helium in our Universe these extra photons are not as impor- 
tant for the total recombination radiation as those from hydrogen. 
Still for accurate computations of a spectral template one should 
also take these into account, but in this case also aspects related 
to det ails in the L yman a radiative transfer dChluba & Sunvaevl 
l2009d f ah. feedback Jchhiba & Sunvaevll2009bl). and electron scat- 
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tering iRubifio- Martin et alj|2008l : IChluba & Sunvaevll2008al) will 
become important. 

In addition, one should mention, that, for example, in the case 
of decaying particles it is possible that during helium recombina- 
tion much more energy is released than during hydrogen recom- 
bination. In that situation the total emission from helium could be 
increased many times, although the changes to the hydrogen recom- 
bination spectrum may still be small. Similarly, the modifications to 
the ionization history might still be not important for computations 
of the CMB power spectra, so that only the cosmological recombi- 
nation spectrum will allow us to put constraints on possible extra 
energy release in this case (see Sect.|4]for more discussion). 



3.3 Dependence of the cosmological recombination spectrum 
on ionization and excitation efficiencies g* 

In the previous sections we have assumed that all the are 
given by the expression l|7}, which was initially suggested by 
IChen & Kamionkowskil d2004l) . In this section we demonstrate that 
the DM-induced contribution to the cosmological recombination 
spectrum is sensitive to the branching of the deposited energy into 
heating, ionization and excitations. This emphasizes how important 
it is to refine the modelling of the effect of DM annihilations on the 
recombination process, especially when aiming at computing de- 
tailed templates for the recombination spectrum or accurately ac- 
counting for this process in connection with the CMB power spec- 
tra. 



3.3.1 Changes in the recombination spectrum when excitations 
by DM annihilations are not included 

As a first case, we will assume that excitations caused by DM an- 
nihilations are negligible. If additional excitations are mediated by 
photons, such approximation will be more appropriate. It is clear 
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Figure 6. Contribution to the total H i cosmological recombination spec- 
trum which is induced by dark matter annihilation for /jeo = 2 X 10~ 22 and 
different choices of g". We included 20 shells for hydrogen and helium in 
our computations. Like in Fig. [4] we give the absolute difference, A/,j m (v), 
in the total, pres ent-day Hi recombination spectrum with respect to our 
reference model jRubino-Martfn et al. 2008) which does not include DM 
annihilations. For the dashed-dotted line we used g? as given by Eq. jTJ, 
for both hydrogen and helium. In the case 'no excitations' we excluded the 
DM-induced excitations of hydrogen and helium. For the solid curve we 
use the expression (8j for all g". 
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Figure 7. Ionization efficiency, g^, for the standard ionization his- 
tory. The solid like corresp onds to the approximation proposed by 
IChen & Kamionkowskil 120041), while the dashed curve is motivated by the 
work of lShull & van Steenberd Jl985h . We also marked the region around 
the last scattering surface (z ~ 1 100), and where most of the recombination 
lines from hydrogen (z. ~ 1400) are produced. 



that this will reduce the amount of extra low frequency emission, 
since due to the strong coupling of e lectrons in the 2p state to the 
3d and 3s state, and the continuum JChluba & Sunvaevll2009e] |cT). 
every additional excitation to the second shell also leads to some 
additional ionizations and transitions among highly excited states. 

In Fig. [6] we show the DM-induced contribution to the total 
H i cosmological recombination spectrum when neglecting excita- 
tions by DM annihilations. At low frequencies the overall ampli- 
tude of the distortion is reduced by ~ 15% - 20%, while at high fre- 
quencies the distortions are about two times lower than in the case 
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which includes excitations. Also the changes close to the B aimer 
and Paschen features reduced by a factor of ~ 1.5. For the same 
/dfo one therefore finds a smaller admixture of the DM-induced 
signal to the normal cosmological recombination spectrum. This 
implies that the changes in the width and position of the recombi- 
nation lines will be smaller when excitations are not efficient. 

It is clear that part of the difference can be compensated by 
increasing the effective value f^eg, but since also the relative ampli- 
tudes of the features are affected (e.g. low to high frequency con- 
tributions), a differential signal remains. This in principle should 
allow to determine how efficient excitations from DM annihilations 
are. As mentioned in Sect. 12. 3.T1 this will depend on details of the 
DM model and the annihilation channels that are important. 

3.3.2 Direct dependence on the redshift scaling ofg* 

As pointed out bv lChen & KamionkowskH d2004l) . and as also men- 
tioned in Sect. |2.3.T1 the approximation Q is very crude and does 
not capture most of the real dependencies of the branching of the 
deposited energy into heating, excitations and ionizations on the 
ionization degree, density of the plasma, helium abundance and ex- 
pansion rate of the Universe. However, the modifications to the re- 
combination spectrum do depend on the scaling of g" with redshift. 
As another example, we therefore also compute the changes in the 
recombination spectrum using the approximation instead of Q. 

In Fig. [7] we present the comparison of these two approxi- 
mations for the stan dard ionization history computed with Recfast 
dSeager et al]| 19991) . As one can see, the ionization efficiency rises 
much slower when usi ng the expression ([8), which i s more closely 
based on the work of IShull & van Steenberd dl985l) than Eq. Q. 
At z ~ 1 100, i.e close to the maximum of the Thomson visibility 
function, g"f hu " is practically two times smaller than g^ he °- This 
implies that ionizations will become important significantly later, 
implying that also the maxima of the DM-induced spectral features 
should be shifted towards higher frequencies (see Sect. 13.2.11 for 
explanation). Furthermore one expects that the overall amplitude 
of the additional distortions should be smaller, since the energy de- 
position rate dBa/ dt decreases with redshift. 

In Fig. [6] we present the DM-induced distortion when using 
-h, u j or a jj ~ a j n( j ee( j me maxima of the distortions are shifted 

° ion ° i 

by Av/v ~ 10% - 20% and the overall amplitude reduced by 
20% - 50% at different frequencies. This shows that the cosmolog- 
ical recombination spectrum is rather sensitive to the detailed time- 
dependence of the ionization and excitation efficiencies. Given the 
large uncertainty in these functions in the current computations, 
it will be important to refine the modelling of the energy depo- 
sition by DM annihilations in this respect. This will also be im- 
portant in connection with precise computations of the ionization 
history and CMB power spectra, and the obtained limits on models 
of annihilating DM using current CM B data (e.g. lGalli et alj|2009l ; 
ISlatver et alj2009l ; lHuetsi et alj|2009l) . 

3.4 Effect of DM annihilations including the 
time-dependence of / d 

In the previous sections we have seen that the shape of the addi- 
tional spectral distortion did not change very much when increas- 
ing the DM annihilation efficiency (see Fig. [4j. However, neglect- 
ing the effect of excitations or changing the redshift dependence 
of gf did lead to some notable modifications in the shape of the 
distortions (cf. Fig.|6). Similarly, one expects that including the ad- 
ditional time-dependence in the energy deposition rate will affect 



60 r— 

30 - 

10 - 



20 shell Hydrogen & 20 shell Helium atom 

5; : — ' — ~i 1 1 1 1 — 

: reference model 

• « ,\\ : ( 1 f (1 =2x 10' 17 ,r&K2(HU 

"\\ ^f = 2x 10"' ? .C&K 2004, no excitation 

\: \ f A t Q = 2 x 10 22 , Shull 1985 

7 ^ 

End of Hal 




200 400 600 



1000 1200 1400 1600 1800 2000 2200 2400 2600 2800 



Figure 8. Effect of dark matter annihilation on the free electron fraction 
for /d^o = 2 X 10~ 22 and different choices of g?. We included 20 shells for 
hydrogen and helium in our computations. The curves in the upper group 
show the relative dif ference in the free electron fraction in comparison to 
the reference model ( Rubino-Martfn et al .120081) . The shaded area indicates 
the region around the Thomson visibility function, which defines the last 
scattering surface. For the dashed-dotted line we used gf as given byEq. 0, 
for both hydrogen and helium. In the case 'no excitations' we excluded the 
DM-induced excitations of hydrogen and helium. For the solid curve we 
use the expression {8j for all g". 



the distortions. Here we want to compare the distortions for dif- 
ferent models of annihilating DM, however, as we show below the 
model-dependence of /d intro duces only rat h er sm all differences. 

Inspecting the results of ISlatver et al.l J2009T) for the energy 
deposition efficiencies, /d(z), for different models of annihilating 
DM, one can see that in most cases only the overall amplitude of 
/ d is changing, while the shape is very similar, resembling the one 



for XX 



e + e DM annihilation with M y 



lOOGeV. This leads 



to a strong degeneracy, since such changes in the overall ampli- 
tude can be compensated when allowing for appropriate (constant) 
boost factors to the annihilation cross-section, e.g. motivated by the 
effect of Sommerfeld enhancement. Therefore, we expect that in all 
these cases the shape of the DM-induced spectral distortion will be 
very similar, and that only the amplitude will depend on the specific 
model via an overall efficiency factor. 

The largest differences in the shape of /d(z) can be found 
for DM models that annihilate via the channels XX ~ * e+e ~ with 
M x = 1 GeV and M x - lOOGe V, and yx ~» /" + /"~ D M annihilation 
with M x = 1 GeV (see Fig. 4 in lSlatver etal]|2009l) . The functions 
fi which we used to represent these cases are shown in Fig. [9] The 
curves were computed applyin g the fitting formul ae for the dif- 
ferent redshift ranges given bv lsiatveret~aLld2009l) and smoothly 
connecting them at z ~ 170, and linearly between z ~ 1470 and 
z ~ 2500. One can see, that around the time of maximal photon 
production by helium (z ~ 2200), all of these functions are more or 
less constant at a level of ~ 1%. Similarly, during the time of photon 
release by hydrogen (z ~ 1300 - 1400) all the functions / d are only 
weakly dependent on time. The largest time-dependence is seen for 
the case^p — » e + e~ with M x = 100 GeV. In all shown cases, the 
strongest variations appear at much lower redshifts (z < 600), so 
that from the differences in /d one does not expect any important 
changes to the shape of the CMB distortions induced by the dif- 
ferent DM annihilation models. The dependence on the branching 
into heating, ionizations and excitations is much more crucial. 
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Figure 9. Energy deposition efficiency, / d , for different models of annihi- 
lating DM. We also marked the region around the last scattering surface 
(z ~ 1100), and where most of the recombination lines from hydrogen 
(z ~ 1400) and neutral helium (z ~ 2100) are produced. 



by the value at high redshift, fixing f£ m 6o = 2 x 10" 22 . One can 
see that the results in the cases of XX ~~ * e+e an( ^XX ~ * A'V DM 
annihilation with M x = 1 GeV are extremely similar. This is simply 
a consequence of the tiny differences in the redshift dependence of 
/d(z)// d Jm for the two models at z ~ 1200 - 1300 (cf. Fig.|9}, from 
where most of the extra H i distortions are coming from. The small 
difference with respect to the case f&(z) = const is due to the fact 
that at z ~ 1300 - 1400 in both cases the value of f d (z)/f^ m < 1. 
This also explains the additional difference in the overall amplitude 
seen for the model with^ — » e + e~ with M x = 100 GeV. 

If in addition we rescale the distortions so that they all coin- 
cide!] at v ~ 200 GHz, then we find that percent-level differences in 
all cases remain. This shows that the model-dependence in connec- 
tion with /d(z) is very small and hence the DM-induced distortions 
will not be very sensitive to the total mass and magnitude of the 
annihilation cross-section by means of this quantity. This is also 
true f or the effects on the CMB power spectra (e.g. ISlatver et al.l 
120091) . however, we expect that more detailed computations of the 
ionization and excitation efficiencies could reveal more pronounced 
dependencies on the specifics of the DM annihilation model. 
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Figure 10. Contribution to the total H i cosmological recombination spec- 
trum which is induced by dark matter annihilation for f^ m £o = 2 X 10~ 22 
and different choices of fi(z)/ f^ m - We included 20 shells for hydrogen and 
helium in our computations, with g d given by Eq. (7). Like in Fig. [3] we 
give the absolute difference, A/d m (v), in the total, pre sent-day H i recombi- 
nation spectrum with respect to our reference model jRubino-Mart m et alj 
|2008|) which does not include DM annihilations. For the dashed-dotted line 
we used fi(z)lf^ m = const, while for the others we used the different fi(z) 
shown in Fig. [9] 

The same statement applies to DM-induced modifications in 
the ionization history and the CMB pow er spectra, and was al- 
ready pointed out by ISlatver et al.1 J2009t) . Only the overall am- 
plitude and its conversion into DM mass and annihilation cross- 
section is affected. However, since of the CMB anisotropies f A (z) at 
z ~ 1 100 matters, while for the Hi spectral distortions the value at 
z ~ 1200- 1400 is more important, a combination of both the CMB 
energy spectrum and power spectra could help shedding additional 
light on the time-dependence of /d(z). Still, in the considered case 
the variation of /d(z) during this time is rather small. 

In Fig. [TO] we finally show the changes in the cosmological re- 
combination spectrum for different models of annihilating DM. To 
make the results comparable in amplitude, we have rescaled f&(z) 



4 DISCUSSION 

In the previous sections we have focused on illustrating the main 
dependencies and effects in connection with the changes that are in- 
troduced by DM annihilations during cosmological recombination. 
In this section we want to go slightly beyond a purely theoretical 
study and include current limits on the possible DM annihilation 
rate into our considerations. We also wish to extend the discussion 
to more general cases of energy release and explain why the cosmo- 
logical recombination spectrum could allow us to learn something 
about non-standard thermal histories in cases to which the ioniza- 
tion history and CMB power spectra are not sensitive. 

4.1 CMB power spectra versus the cosmological 
recombination spectrum 

For models which still seem to be all owed by current CMB 
data from Wmap dKomatsu et~aU l2009h. at 95% c.l. one has 
fd <o-v> ^¥ S 3.6x 10~ 25 cm 3 s" 1 jGalli et al.l2009llSlatver et~ai1 
120091) . In our parametrization, according to Eqs. ((2) and {3), this 

translates into f^eo < 2 x 10~ 23 [t]t] ■ F° r m ' s value the changes 
in the cosmological recombination spectrum are expected to be of 
the order of percent only (cf. Fig.[3jand Fig. [4}. Unfortunately, this 
would make it rather hard to learn something in addition about DM 
annihilations from the cosmological recombination spectrum. So 
why should one try to measure the cosmological recombination 
spectrum when the constraints obtained with the CMB temperature 
and polarization power spectra are already so strong? 

First of all, one should emphasize that for the allowed values 
of fi (crv) also the changes in the CMB temperature and polariza- 
tion power spectra are fairly small: For / d (crv) = 2 x 10~ 23 they 
reach ~ 4% - 5% at high multipoles (see Fig. II It. where the main 
trend is completely featureless, and the variable part has an ampli- 
tude of only ~ 1%. In this context it is particularly important that 
due to cosmic variance percent-level correction to the CMB power 
spectra only become statistically significant at high multipoles. 

2 We arbitrarily chose this frequency. 
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Figure 11. DM-induced changes to the CMB temperature and polarization 
power spectra (upper panel) and the H i cosmological recombination spec- 
trum (lower panel) for f\' m £o = 2 X 1 0~ 23 . We used Cmbeasy to compute 
the changes in the CMB power spectra jDoranli2005l) . 



Secondly, it is clear that the DM-induced modifications in the 
CMB power spectra are strongly degenerate with the values of the 
scalar spectral index n s , the baryon density Q. b , and o ~g (e.g. see 
Padmanabhan & Finkbeinen 120051 : [slatver et all [20091 : iGalli et~ai] 
20091). Here also previously neglected ph ysical processes (see 
Fendt et al. 20091 ISunvaev & Chluba 2009, for overview) are im- 
portant, since for Pl anck they can lead to signifi cant biases to the 
values of n s and Q b jRubino-Martm et"al or otherwise an- 

other confusion to the possible signatures of DM annihilations. 

Furthermore, at high multipoles the signal s related to SZ- 
clusters (e.g. see iMolnar & Birkinshawl |2000|: iRefregier et al .1 
l200d:lHolder & Carlstromll200lLlKomatsu & Seliakll2002l) and eiT 
perimental systemat ics (e.g. due to calibrations of the beam, 
IColombo et al]|2009h become very important. 

Disentangling all these components is a mayor problem, and 
this is where the cosmological recombination spectrum could help 
in addition: just looking at the amplitude of the DM-induced 
changes the cosmological recombination spectrum, it is clear that 
in this respect the recombination spectrum in principle is similarly 
sensitive as the CMB power spectra (see Fig. [TTT l. Also, the cos- 
mological recombination spectrum does not depend on the value 



of the scalar spectral index n 8 , so that such confusion is already 
excluded. In addition, the cosmological recombination spectrum is 
not limited by cosmic variance: no statistical comparison of the 
measured energy spectrum with some ensemble of Universes is in- 
volved. One would investigate the recombination spectrum for our 
particular realization of the Universe, only encountering small fluc- 
tuations of the cosmological parameters in different directions of 
the sky. It even is possible to take these corrections into account 
in the computations of the cosmological recombinations spectrum, 
but for standard cosmological model the effects are expected to in- 
troduce changes that are smaller than A/ v // v ~ 10~ 4 - 10" 3 . 

To observe the cosmological r ecombination spectrum no 
absolute measurement is necessary dChluba & Sunvaevl l2008bl: 
ISunvaev & Chlubal l2009h . Due to its very peculiar frequency- 
dependence, it should therefore in principle be possible to separate 
the cosmological recombination spectrum from foregrounds and in- 
strumental signals. One can use particularly clean patches on the 
sky, observing with wide-angle horns. Also the cosmological re- 
combination signal should be polarized at a very small level only, 
providing another possibility to discriminate it from other signals. 
By measuring the positions, and width of several spectral features 
in the cosmological recombination spectrum with high precision 
one should therefore in principle be able to extract valuable in- 
formation on possible sources of extra ionizations and excitation 
during recombination. 



4.2 Why the cosmological recombination spectrum could 
teach us something in addition 

Until now we have only considered the possibility of continuous, 
very extended energy release caused by DM annihilations. How- 
ever, as an example, for energy injection due to long-lived decaying 
particles additional aspects become important: First of all, in these 
cases most of the energy will be released over a characteristic time 
Az/z ~ 30% - 40%. This implies that for decaying particles the 
possible changes to the ionization history or cosmological recom- 
bination spectrum could be significantly more sharp. Secondly, if 
the life-time of the particle is shorter than ~ 380000 years, then 
most of the energy will be released before the maximum of the 
Thomson visibility function. In that case the effect on the CMB 
temperature and polarization power spectra could be much smaller 
or even negligible, while the changes in the cosmological recombi- 
nation spectrum would still be very strong. 

For example, if energy were released at ~ 1 8 000 years after 
the big bang (corresponding to the time of He in — > He ii recom- 
bination) then the CMB power spectra would practically remain 
unchanged, while at the same time due to extra ionizations and ex- 
citations the number of photons produced by Heii could increase 
significantl)[|. Similarly, if the energy is released between the two 
epoch of helium recombination or during the recombination of neu- 
tral helium, the CMB power spectra will mostly be unaltered, while 
a lot of extra emission by He ii could be induced. The cosmological 



3 Such behavior was seen in connection with the possibility of pre- 
recombinational energy release, where the contributions from Heii to 
the recombination radiatio n became as large as those from hydrogen 
IChluba & Sunvaevll2009"3) . However, we note that in the case considered 
here, there is no significant >'-type distortion produced by the energy injec- 
tion, but as mentioned in the introduction, one can still have a very impor- 
tant effect on the cosmological recombination spectrum, since we allow for 
direct ionizations and excitations of neutral atoms. 
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recombination spectrum therefore may allow us to check for ex- 
tra sources of ionizations or excitations at times where the CMB 
temperature and polarization power spectra are not sensitive. 

It is also important to mention, that any other modification 
of the recombination process, e.g. caused by the variation of fun- 
damental constants like the fine-structure constant a, or Newtons 
gravitational constant G (see lScoccola et alT2 009; Gall i et al.l2009l 
for references), will also lead to changes in the shape and positions 
of features in the cosmological recombination spectrum. Observing 
such changes may therefore provide another very clean and direct 
way to test non-standard physics at early stages of the Universe. 



4.3 Towards detailed templates for the cosmological 
recombination spectrum 

There is no principle difficulty in computing the cosmologi- 
cal recombination spectrum with ~ 0.1% accuracy, also includ- 
ing the effects of possible energy or particle release in the re- 
combination epoch, even for more general cases. However, in 
this case, one would ha ve to take other p r eviously neglected 
physi cal processes (see iFendt et alj 120091 ; ISunvaev & Chlubj 
I2009L for overview) into account in addition. Here in par- 
ticul ar those connected with the Lyman a rad i ative trans- 
fer dChluba & Sunvaevl l2009ellal: iHirata & Forbes! l2009h and 



two-photon transitions jghluba & Sunvaevl l2008d : Ifflratal 120081 ; 
IChluba & Sunvaevl 1200937 will be very important, because they 



10% 



(see comments in IChluba & Sunvaevl 


2009d). Also the effect of 


electron scattering will be important ( 


Rubino-Martm et alj|2008fc 



helium. 

Furthermore, the computations of the heating, ionization and 
excitation efficiencies will also have to be refined in order to 
study detailed model-dependencies (for additional comments see 
Sect. [23}, and as we have seen in Sect. 13. 31 the difference in the re- 
combination spectrum can be large. Also one should include more 
shells for hydrogen and helium into the computations, since the 
total amplitude of the additional distortions is still expected to in- 
crease, especially at low frequencies. 

In addition, the emission due to He n should be taken into ac- 
count. Given that there is a very extended period between the re- 
combination of He 11 at z ~ 6000 and He i at z ~ 2200, the total 
amount of He n emission induced by DM annihilations is expected 
to be very important, possibly even exceeding that from neutral he- 
lium. Such calculations are beyond the scope of this paper, but we 
plan to investigate these aspects in more detail in the future. 

However, we note that currently it will probably be more im- 
portant to investigate more general observational prospects in con- 
nection with a measurement of the cosmological recombination 
spectrum, including possible foregrounds and systematics. Also, it 
will be very important to understand, which frequency bands will 
be most useful and sensitive to changes in the cosmological pa- 
rameters or energy injection. It is obvious that measuring the cos- 
mological recombination spectrum will be very challenging, but on 
the other hand there could be a lot to learn from this. In particular, a 
combination of the CMB temperature and polarization anisotropics 
with the CMB energy spectrum could open a way to further tighten 
CMB-based constraints on standard and non-standards aspect of 
our cosmological model. 



5 CONCLUSIONS 

Assuming that DM is annihilating throughout the history of our 
Universe, we have demonstrated that the cosmological recombina- 
tion spectrum in principle is sensitive to the branching of the de- 
posited energy into heating, ionizations and excitations. If energy 
only goes into heating of the medium (without leading to some sig- 
nificant primordial ft or y-type CMB distortion), the recombination 
spectrum in practically not affected, while extra ionizations and ex- 
citations lead to modifications in both the contributions from hy- 
drogen (cf. Fig. [3} and helium (cf. Fig.[5}- 

We have shown, that the overall amplitude of the DM-induced 
spectral distortions depends on the total amount of ionizations 
and excitations (Fig. [3J at z ~ 1200 - 1300 for hydrogen, and 
z ~ 2000-2400 for helium. Furthermore, the relative importance of 
DM-induced excitations and ionizations, and the time-dependence 
of their efficiencies determine the exact shape and position of the 
additional CMB spectral distortions from the recombination epoch 
(see Fig. [6). Since these efficiencies depend on the considered 
model for the annihilating DM, or more generally, the process that 
produced the additional ionizations and excitations (e.g. decaying 
particles), by measuring the cosmological recombination spectrum 
in several spectral bands in principle one should be able to place 
additional constraints on possible energy release during cosmolog- 
ical recombination. Given the rather strong dependence of both the 
changes to the cosmological recombination spectrum and the CMB 
temperature and polarization power spectra on these efficiencies, 
for precise predictions in connection with possible extra release 
of energy during recombination it will be important to refine the 
computations of these efficiencies in the cosmological context (see 
Sect. |2.3| for more comments). 

Although for currently allowed values of the effective DM an- 
nihilation rate the changes to the cosmological recombination spec- 
trum are of the order of percent, we have argued that there are sev- 
eral reasons to believe that one could learn something in addition 
by studying the signals from the recombination era (see Sect.[4]l. In 
particular, the cosmological recombination spectrum is expected to 
be sensitive to cases of energy or particle release (e.g. before the 
maximum of the Thomson visibility function), by which the CMB 
temperature and polarization power spectra are not affected (see 
Sect. 14. 2t . Beyond a check that recombination has occurred as we 
think it has, observing the cosmological recombination radiation 
would therefore allows us to directly check our understanding of 
standard and non-standard physical processes happening at about 
260000, 130000 and 18 000 years after the big bang. 
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APPENDIX A: ESTIMATES FOR THE PRIMORDIAL 
CMB SPECTRAL DISTORTIONS FOR THE CONSIDERED 
MODELS OF ANNIHILATING DM 

A significant fraction of the total energy released due to DM anni- 
hilation is going into heating of the medium, in particular at red- 
shifts well before the recombination epoch, where f, w /d and 
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/ion ~ /ex ~ 0. It is clear that this will lead to some (tiny) primordial 
distortion of the CMB, since at redshifts z < 2x 10 6 the thermaliza- 
tion process stops being 100% efficient (e.g.|lllarionov & Svunvaevl 



tion process stops being 1UU% efficient (e.g. Ill 
ll975bl : lBurigana et aljl99ll ; lHu & Silkll 19931) 



One can estimate the expected primordial distortion of the 
CMB by computing the total change in the energy density of the 
CMB, Apy/p-y = J — dt, that is caused by the deposition of en- 
ergy in the DM annihilation process. For a p and y-type spectral 
distortion one then expects yL » 1.4Ap r /p r and y « jAp y /p y (e.g. 
Illlarionov & Svunvaev | Il975bl). For s-wav e annihilation ((crv) = 
const) one finds (cf. lMcDonald et alj200lh 



3x 10~ lu /d, 



100 GeV 



0.13 



(crv) 



3 X 10- 26 cm 3 /s 



(Al) 



where /djim is the limiting energy deposition efficiency at large 
z. This is only a rough estimate, but since the final value for 
fi is so far below the c urrent upper limit \ y\ < 9 x 10~ 5 
obtain ed with Cobe/Firas dMather et al.l 1 19941 ; iFixsen & MatheJ 
|2002|). for our purpose it is sufficient. The effective ^-parameter 
i Zeldovich & Sunvaev|[l969h is of the same order. 

Although for y ~ 10~ 7 - 10~ 6 a significant amount of pho- 
tons is gene rated in the pre-recombina tional epochs of hydrogen 
and helium dChluba & Sunyaev! l2009dl) . inducing interesting nar- 
row spectral features in the cosmological recombination radiation, 
at the level of /j or y ~ 10~ 10 this process can be completely ne- 
glected. In our computations we can therefore assume that the am- 
bient CMB radiation field i s given by a blackbody with present-day 
temperature T = 2.725 K dFixsen & Mathedl2002l) . 
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